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ABSTKACT 

The  purpose  of  this  paper  is  to  collect  certain  wave  propagation 
results  in  two-snd-one-half  distensions  —  defined  as  three  diaensional 
propagation  in  a  aediua  that  has  variations  in  two  diaensions  only.  The 
results  of  interest  are  for  sources  and  receivers  in  the  plane  deterained  by 
the  two  directions  of  paraaeter  variation.  The  objective  of  this  work  is  to 
reduce  the  analysis  of  the  in-plane  propagation  to  tvo  diaensional  analysis 
while  retaining  —  at  least  asyaptotically  —  the  proper  three  diaensional 
geoaetrical  spreading.  Ve  do  this  for  the  free  space  Green's  function  and 
for  the  Kirchhoff  approziaate  upward  scattered  field  froa  a  single 
reflector.  In  both  cases,  we  carry  out  a  derivation  under  the  assuaption  of 
a  background  velocity  with  two  diaensional  —  c(x, z)  —  variation/  we 
specialize  the  results  to  a  constant  background  velocity  and  a  depth 
dependent  background  velocity.  For  the  convenience  of  the  user  we  have 
included  a  glossary  and  two  tables  of  equation  nuabers  to  help  in  finding 
specific  results. 


aaplitode  of  ray  theoretic  Green'*  function  in  3D,  constant 
density. 


aaplitode  of  ray  theoretic  Green's  function  in  2D,  constant 
density. 

ray  parasieter  labeling  rays,*  polar  angle  of  initial  direction 
in  (x, y) -plane. 

aaplitode  of  ray  theoretic  Green's  function  in  3D,  variable 
density. 

aaplitode  of  ray  theoretic  Green's  function  in  2D,  variable 
density. 

ray  paraaeter  labeling  rays /  axiauthal  angle  with  respect  to 
zt  becoaes  polar  angle  in  (x, s)  plane  when  y  =  0. 

background  velocity. 

background  velocity  depending  on  x  alone. 

propagation  velocity  below  (above)  reflector,  S. 

Jacobian  of  aapping  via  rays  in  3D,  d(xt,  x,,  xg  )/3(o,  a.  ft) . 

Jacobian  of  aapping  via  rays  in  2D,  3(xt, x, )/8(o,  0),  a  *  0. 

index  of  refraction  with  c(z)  background  velocity, 
n(x)  «  c($,)/c(z) . 

Vt,  p*p  *  1/c*. 

ray  theoretic  reflection  coefficient  defined  by  (46), 
specialized  to  constant  background  velocity  in  (62). 


(x1  -  $4)  +  (x,  -  {,) 


variable  density. 

ray  paraaeter  for  which  i1  *  1/c1  along  a  ray 


travel tiae. 


Vb**^ 


upward  scattered  wave  from  a  reflector  S  due  to  a  point 
source  at  £  evaluated  at  3. 

Cartesian  coordinates. 


(x,  y,  x) 

(*a»xa»**)  Cartesian  coordinates, 

x  Cartesian  coordinates. 

£  (£x>0.£t),  source  point  for  Green’s  function. 

Zp  turning  point;  n(zp)  =  sinf). 
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icnm 


Hie  gathering  of  seiaaic  data  over  a  line  on  the  surface  of  the  earth 
or  ocean  —  rather  than  over  a  planar  or  areal  array  —  is  still  a  pre¬ 
eminent  aethodology.  Consequently,  aethods  of  reflector  imaging  and  earth 
parameter  estimation  —  migration,  structural  inversion,  seismic  inversion 
—  must  be  based  on  models  of  the  substructure  and  vave  propagation  which 
can  produce  meaningful  results  from  this  dimensionally  constrained  data. 


Finite  difference  wave  equation  migration  [Claerbout,  [1976]  deals  with 
this  difficulty  by  assuming  a  two  dimensional  model,  both  for  the 
substructure  and  the  wave  propagation.  The  two  dimensions  are  the  depth  — 
x  or  i,  —  and  the  transverse  range  along  the  line  of  the  seismic 
experiments  x  or  x4.  While  it  might  be  reasonable  to  assume  that  the  out- 
of-plane  variationa  are  email  enough  to  be  neglected,  the  implementation  of 
a  two  dimensional  vave  equations  guarantees  that  the  amplitudes  generated 
will  be  progressively  more  inaccurate  with  traveltiae  due  to  improper 
characterisation  by  the  two  dimensional  wave  equation  of  the  geometrical 
spreading  effects  in  three  dimensions. 


On  the  other  hand,  lirehhoff  migration  [Sohneider,  1978],  k-f  migration 
:?tolt,  1978],  and  Born  Inversion  [Cohen  and  Bleistein,  1979],  among  others, 
pr.'perly  allow  for  three  dimensional  propagation  of  waves  —  within  die 
accuracy  of  the  assmnptions  of  these  aethods.  However,  neither  these 
methods  nor  any  other  can  treat  three  dimensional  variation  in  the 
substructure  when  only  a  linear  surrey  is  carried  out  on  the  surfaee.  The 
limitation  of  the  data  is  accomodated  by  assuming  that  the  earth  parameters 
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▼ary  only  with  the  in-plane  variables  —  (x,  z)  —  and  are  independent  of  the 
ont-of  plane  variable,  y  or  x»> 


Thus,  the  earth  variations  are  essentially  two  dimensional  — 
cylindrical  in  three  dimensions  —  while  the  propagation  is  three 
dimensional.  A  few  years  ago,  I  heard  a  presentation  by  6.  Hohmann  on  an 
electromagnetic  problem  with  the  same  type  of  geometry.  He  called  this  case 
two-and  one-half  dimensional  —  2*7,0.  I  have  adopted  that  terminology,  as 
well . 


In  2 V ,D  one  is  almost  always  interested  in  only  in-plane  values  of  the 
wave  field/  that  is,  values  of  wavefield  with  y  or  x,  equal  to  zero. 
Certainly  one  needs  only  in-plane  values  of  the  earth  parameters.  Thus, 
except  for  the  effect  of  three  dimensional  spreading,  the  problems  of 
interest  are  essentially  two  dimensional. 

In  research  in  our  own  group,  we  are  constantly  in  need  of  in-plane 
evaluations  of  24/,D  asymptotic  wave  fields  both  for  inversion  algorithms 
and  for  direct  modeling.  Ve  find  that  we  can  reduce  such  wave  fields  to  2D 
wave  fields  multiplied  by  appropriate  factors,  also  determined  totally  in 
2D.  The  purpose  of  this  note  is  to  derive  and  tabulate  some  of  those 
results.  The  starting  point  for  these  results  is  either  three  dimensional 
ray  theory  or  a  generalized  version  of  the  Kirchhoff  approximation  for  the 
upward  scattered  wave  for  a  single  reflector.  These  results  are  available 
from  many  sources.  However,  I  shall  use  my  own  booh  [Bleistein,  19841  — 
referred  to  below  as  MNVP  —  because  I  can  find  those  formulas  most  easily 


there 


We  chall  state  results  for  a  c(x. a)  background  velocity  and  specialise 
those  to  constant  background  velocity  where  the  reaults  becoae  explicit  and 
to  a  c(x)  background  velocity  where  the  results  are  soaewhat  aore  explicit 
than  in  the  general  case.  For  aost  of  the  paper,  we  discuss  the  case  of 
constant  density  and  variable  soundspeed.  However,  in  Section  5  we  describe 
the  necessary  aodif ications  for  variable  density.  For  the  convenience  of 
the  user  I  have  included  a  glossary  of  the  notation  of  this  paper  and  two 
tables  of  ^ution  nuabers  to  help  in  finding  specific  results. 


2.  ASYMPTOTIC  2V,D  GKESN'S  FUNCTION 


Ve  consider  the  following  problem  for  n: 


T»+^n-  -6(xt-  t4)  6<x#)  5(xt-  5f) 


In  thie  equation,  V*  is  the  Laplscisn  in  the  variables  x  =  (x,  y,  z)  = 
(x1# X*.  xt)  and  c  =  c(x, z)  =  c(xx,xt).  We  use  the  convention  that  z  =  x,  is 
positive  downward.  We  propose  to  solve  this  problem  by  ray  methods  [Keller, 
1958/  Lewis  snd  Keller,  1964].  Oar  objective,  here,  will  be  to  reduce  this 
ray  theoretic  problem  in  three  dimensions  to  a  problem  in  two  dimensions 
when  x,  =  0  while  still  retaining  the  correct  solution  of  the  three 
dimensional  problem  by  appropriate  adjustment  of  the  equations  and  scaling 
of  the  solution. 


We  assume  that 


u(x,*w)  ~  A(x)  exp{iur(x) } 


with  x  a  solution  of  the  eikonal  equation. 


(Vc) *  =  l/c*(x)  , 


and  A  a  solution  of  the  transport  equation 


2  WA-Vt  +  AV  t  -  0  . 


(See  MNVP,  Section*  8.2  end  8.3,  for  details  of  the  derivation  of  these 


equations  and  results  below  related  to  ray  theory.)  The  functions,  t  and  A, 
nust  also  satisfy  the  conditions 

x  ■  0,  for  x  =  £  =  (Si.0,5,)  (5) 

Jx-£ |  A  -)  l/4n,  as  x  ->  £  .  (6) 


Ve  shall  solve  (3),  (5)  by  the  method  of  characteristics.  The 

characteristic  equations  or  ray  equations  are 


*  -  p  . 


g  *  Vt,  2*  *  l/c*(x) 


1 

P  =  y 


V(l/c*(x))  , 


(7) 


X  =  l/c*(x)  , 


The  solution  we  seek  is  the  conoidal  solution  for  which  the  rays  enanate 
froB  a  point,  £.  We  require  initial  data  at  a  *  0  for  each  of  the  seven 
unknowns  in  (7).  The  data  for  t  is  given  in  (5)/  since  the  rays  eaanate 
froa  £,  the  initial  values  of  x  are  known.  For  the  conoidal  solution,  the 


initial  value*  of  £  are  not  known,  except  that  they  most  be  constrained  by 
the  eikonal  equation,  (3),  itself:  p*  =  1/c*.  Thus,  there  are  two  other 
paraaieters  which  serve  to  label  a  ray  by  its  initial  direction.  Therefore 
the  initial  data  or  ray  data  for  (7)  is 

x(0)  =  J 

j>(0)  =  l/c(£)  (cos  a  sin  0,  sin  a  sin  0  ,  cos  0  )  .  (8) 

t(0)  -  0 

The  solution  of  ordinary  differential  equations  (7)  and  the  ray  data  (8)  is 
a  family  of  rays  with  parameter  a,  distinguished  from  one  another  by  the 
choice  of  the  parameters  a  and  0.  Along  each  ray,  the  value  of  t  is  known, 
as  well. 

The  specific  assumption  that  c  is  independent  of  xa  allows  us  to  obtain 
a  part  of  this  solution  in  closed  form.  Prom  (7),  the  equation  for  p,  is 


(9) 


Thus,  p3  i«  constant  on  each  ray,  that  is,  independent  of  a.  From  the  ray 
data  in  (8),  we  know  the  value  of  that  constant. 


Ps 


sin  a 


sin  0  . 


(10) 


With  p*  independent  of  o,  the  equation  for  xs  in  (7)  now  becomes 


particularly  staple.  Taking  account  of  the  initial  data  in  (8),  we  find 
that 


X*  **  «ina  sin  p  .  (11) 

Ve  remark  that  xa  =  0  for  any  choice  of  (x^  xt)  only  if  we  choose  those 
rays  for  which  sin  a  =  0. 

The  transport  equation,  (4),  can  also  be  written  as  an  ordinary 
differential  equation  in  the  ray  paraaeter,  a.  The  solution  of  that 
equation  is  [MMWP  (8.3.12)]: 


A  = 


4ji 


c(£)J(<r.a.p) 


(12) 


In  this  equation,  J  is  the  Jacobian  of  the  transforaation  froa  x  to  (o,a,p) 
via  the  solution  of  the  ray  equations,  (7)  and  (8): 


J  = 


a(W,s) 
d(a,  a,  p) 


(13) 


Hie  foil  owing  changes  in  notation  from  the  refernce  have  been  aade  here: 
(i)  Yi  ->  fi.  (ii)  p(x#)  ->  1/c  (£);  (iii)  p(x(ff))  J(o)  ->J(o,a,p).  For  the 
last  of  these,  we  have  need  (8.2.29)  with  X  =  1  to  redefine  J  as  in  (13). 

The  solution  (2),  with  r  a  solntion  of  the  system  (7)  and  (8)  and  A 
given  by  (12),  is  valid  np  to  the  first  zero  of  J  for  a  M. 


3.  Asymptotic  In-Plane  2V,D  Gim't  Fn«ti 


Ve  shall  now  consider  the  results  of  the  previous  section  for  the  esse 
x,  >  0.  Ve  shall  refer  to  this  case  as  the  in-plane  solution.  As  noted 
above,  setting  x,  ■  0  is  equivalent  to  setting  a  *  0  in  (8)  and  (11). 

Let  ns  consider  the  renaining  equations  in  (7)  and  under  the  assumption 
that  o  »  0: 


xi  "  Pi  * 

x»  m  P.  * 

•  i  a 


*  2  ax. 


Lc  (*t.  *,) 


j  -  i-i_  [ _ l - .1 

2  «>*,  i  *, r 

*  Lc  (xl.  XjJJ 


a1(0)  =  5X  l 

xi  (0)  -  5,  / 


>.<0,-4t£  - 


/  «os  P 

».<0>  'T7T f  ' 


(14) 


C  (x^,  X,) 


x(0)  -  0  . 


This  is  a  closed  system  of  equations  for  which  pj  +  pj  »  l/c*(x),  since 
p,  =  0.  In  fact,  it  is  just  the  system  of  ray  equations  for  solving  the 
eikonal  equation  in  two  dimensions.  Hence,  an  algorithm  for  solving  ray 
equations  in  2D  may  be  applied  to  determine  the  in-plane  rays  and  phase.  I 
do  not  believe  that  this  would  surprise  anyone  who  has  read  this  far. 


Let  u  now  consider  the  solution  (12)  for  the  amplitude  where  a  =  0. 
In  particular,  from  (11),  we  csn  calculate  the  elements  of  the  second  row  of 


the  matrix  arising  in  J,  (13): 


3x 

a 

1 

To  ~ 

c (ff 

3x 

a 

a 

TT  " 

e(l) 

8*a 

a 

30 

ciV 

sin  a  sin 0»  0,  a  *  0  t 

cos  a  sin  0  =  ^  *  °  =  ®  1 

sin  a  cos  0  =  0  .  a  =  0  . 


(15) 


Ve  see  that  only  the  second  element  is  nonzero.  Thus,  the  3x3  determinant 
in  (13)  can  be  reduced  to  a  2x2  determinant. 


J 


sin0  I  ,* 


(16) 


K  - 


d(xJL.xl) 

3(o,0) 


0*0 


(17) 


The  Jacobian,  K,  can  be  computed  directly  from  the  system  (14)  for  the 
solution  of  eikonal  equation  in  two  dimensions.  Ve  use  (16)  in  (12)  and 
obtain 
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A 


1 


US) 


4h 


By  acini  (14)  cad  (18),  then,  ve  obtain  the  in-plane  21/tD  aayaptotic 
Green's  function  totally  in  tens  of  solutions  of  the  2D  ray  equations.  Ve 
resiark  that  it  was  the  simple  fon  of  this  result  which  motivated  the 
particular  choice  of  ray  parameter  a,  for  which  ±*  ■  l/c*(x),  rather  than 
such  alternative  ray  parameters  as  arc  length  on  a  ray  —  for  which  i*  ■  1 
—  or  traveltime  —  for  which  i*  =  c*(x).  We  note,  however,  that  a  is  an 
unphysical  variable  with  dimensions  (length) */ time. 


The  in-phase  representation  of  the  solution  (2)  now  takes  the  fon. 


exp  (iwc) 

u  ~  .  , 

4n  JwT 


(19) 


with  t,  a  and  K  totally  deteninnble  by  in-plane  calculations.  Nonetheless, 
this  solution  has  the  proper  geometrical  spreading  for  three  dimensional 
propagation. 


The  solution  (18)  can  be  more  precisely  related  to  the  solution  of  tile 
two  dimensional  transport  equation.  (The  reader  should  note,  however,  that 
the  latter  amplitude  is  the  response  to  a  line  source  in  three  dimensions.)  -® 

Let  us  denote  that  solution  in  2D  by  As.  The  result  is  given  by  IMVP  ; 

(8.3.30),  but  in  different  notation.  The  correct  result  is 


-V-Vi 


exp  (la/4  »|a  a) 


(20) 


A*  m 


In  tkli  equation,  we  have  also  included  the  proper  power  of  u  at  given  froa 
the  diacnseion  above  (8.3.30)  and  we  have  accounted  for  negative  e,  through 
the  ahaolnte  values  and  the  factor  sgnw. 


By  comparing  (18)  and  (20).  we  find  that 


A 


A» 


exp  {-in/4  sgn »)  . 


(21) 


thus,  if  one  has  available  a  2D  ray  equation  routine  which  calculates  A, 
directly,  one  can  use  this  routine  to  generate  the  in-plane  2X/*D  aaplitude 
by  employing  the  coaplex  scaling  indicated  by  (21). 


I  further  take  the  position  that  this  aultiplier  is  an  adequate  "quick 
and  dirty”  correction  to  apply  to  an  exact  2D  equation  solver.  Hie  reason 
is  that  in  geophysical  applications,  we  are  alaost  always  interested  in  the 
solution  aany  —  at  least  three  —  wavelengths  froa  the  source  in  a  slowly 
varying  background  aediua.  In  that  liait  (high  frequency),  an  asyaptotic 
adjnstaent  will  suffice. 


It  should  also  be  noted  that  all  refraction  and  reflection  effects  of 
the  in-plane  2X/,D  propagation  will  reaain  in-plane.  Thus,  the  sane 
aultiplier  can  be  used  for  a  discontinuous  c(x, z).  However,  one  aust  take 
care  to  aeasnre  o  froa  the  source  point  and  not  froa  the  discontinuity 


surface.  Furthermore,  one  must  then  use  sa  amplitude  A.  which  properly 
aceouats  for  ia-plaae  traasaissioa  effects. 

For  two  special  choices  of  c(x<x)  we  caa  be  more  explicit  about  the 
relatioaship  between  e  aad  the  Cartesiaa  coordiaates.  Ia  the  simplest  case, 
c  *  coastaat  -  c#»  the  ray  equatioas,  (14).  with  the  iadicated  ray  data  caa 
be  solved  explicitly.  First  we  observe  that  ■  J,  -  0  aad  (pt»pt)  is  a 
coastaat  vector  given  by  the  iaitial  data  with  e(£)  -  c, .  We  thea  solve  for 
x4  and  x,  in  (14): 

Xj  ■  {,  +  c,  *o  sinf),  x,  “  {|+  c  cosft.  t  ■  a/e  (22) 

We  now  solve  for  a  ia  the  first  two  equations 

o  *  cfr,  r  -  ^(xJ-?i)*+(xJ-l#)*  ,  (23) 

and  find,  upon  substitution  into  the  third  equation. 

t  *  r/c#  (24) 

Furthermore,  we  can  find  I  by  applying  the  defiaitioa  (17)  to  the  equatioas 
for  the  rays  ia  (22) .  The  result  is 


fe  as a  (23)  sad  (25)  ia  (18)  to  obtain 

A  -  (4»r)-1  .  (26) 

By  inserting  (24)  sad  (26)  ia  (19)  we  obtsia  the  exact  three  diaeasioaal 
Greea's  faactioa,  which  iadeed.  is  eqoal  to  its  leadiag  order  asymptotic 
expaasioa.  This  serves  as  a  simple  check  on  the  method. 

As  a  secoad  example,  let  ns  suppose  that  e  is  a  function  of  z  above, 
c  -  c(z)  ■  c(xt).  In  this  case,  it  is  easier  to  recast  (14)  as  a  system  of 
equations  ia  z  and  express  a  as  a  function  of  z,  as  well.  We  rewrite  (14)  as 


From  the  second  eq net  ion,  here  ve  see  that  px  is  s  constant  on  each 
ray,  given  by  its  initial  value: 


(28) 


This  is  Snell's  lav  for  a  Medina  with  depth-dependent  velocity.  Since 
(Pi,p()  satisfy  the  eihonal  equation  in  the  fora  in  (7),  ve  conclude  that 


P. 


1 

(z) 


- -  sin*0  . 

c  ($,) 


(29) 


We  have  chosen  the  positive  square  root,  here,  because  ve  are  interested  in 
rays  vhich  are  directed  dovmrards  that  is,  t  aust  increase  vith  a.  At  z*Ct, 
the  square  root  is  real  for  all  8.  For  01  <  w/2,  p,  vill  reaain  real  in 
soae  neighborhood  of  We  continue  our  analysis  in  that  region. 


We  can  now  proceed  to  solve  the  reaaining  equations  in  (27) 


It  is  this  lsst  fsctor  which  we  use  in  (18)  or  (17)  or  (21)  to  adjust  the 
in-plane  Green's  function  for  geometrical  spreading  out-of-plane. 

For  computational  purposes/  it  is  desirable  to  express  these  results  in 
terms  of  an  index  of  refraction: 


n(z) 


c(z) 


(31) 


Those  equations  sre 


n(z')dz 


To  complete  the  computations  for  the  Green's  function  for  this  case,  we 
must  determine  K,  defined  by  (17).  To  simplify  the  computation,  we  note 
that 


0(Xz.X.) 

9(o,0) 

9(zi,z)  9(z,0) 


Idx^  9z 


d(p.z)  TToTpT  [afT 


by  application  of  the  chain  rnle  for  Jacobi  determinants.  The  derivatives 
in  the  last  expression  are  partial  derivatives  because  xx,z  are  determined 
in  the  plane  as  functions  of  the  two  variable  0,  z.  We  calculate  9xx/90  from 
the  first  equation  in  (31)#  9z/9 a  is  determined  from  the  last  equation  in 
(28).  The  result  is  that 


r 


■  'i-  '%.■ 1 *<*•'•* 


i 


ra 


9z. 


dz 


3p  *  do 


z 

(«*<,)  -  iSi-t  f  - pl— ' 1  dl  ’ - _ 

W  J  (.*(,-)  -  .i.**)*" 


(34) 


Ve  can  now  use  (18).  (30)  and  (32)  in  (19)  to  write  down  the  asymptotic 
2*7,0  solution  evaluated  in-plane.  The  result  is  given  parametrically  with 
parameter  p.  For  each  coordinate  pair,  (zt.z),  we  determine  p  from  the 
first  equation  in  (32).  Ve  substitute  this  solution  into  the  second  and 
third  equations  of  (32)  to  determine  x  and  a,  and  we  substitute  into  (34)  to 
determine  K.  When  all  of  these  results  are  substituted  into  (19),  u  is 
determined. 

We  remark  that  the  function  K  in  (34)  is  never  equal  to  zero  for 
*  ^  C| .  That  is,  the  downward  propagating  part  of  the  Green's  function  in  a 
c(z)  medium  has  no  caustics.  The  only  pathology  that  could  occur  is  that 
for  some  p,  the  ray  propagates  to  a  depth,  Zp,  such  that  n(zp)  =  sinp  .  If 

we  assume  that  n'(zp)  £  0,  then  K  has  a  finite  nonzero  limit  at  zp  even 

though  the  integral  in  (34)  diverges  for  z  =  Zp.  To  see  why  this  is  so. 

multiply  and  divide  by  n'(z')  under  the  integral  sign  and  integrate  by 

parts.  The  result  is 


J^kai  -  m.  — *  *-  * 
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*  mm 


^|n*(z)  - 


sinf) 


cot  P 


n(z) 


«(?,) 


n'(z)  i|n*(z)-sin*p  n'({,)  >-sin*0 

(35) 


5, 


n( z  *)  1  _ dzj _ 

t|n*  (z  *)  -  sin*p) 


The  integral,  here,  converges  even  when  z  =  Zp.  Thnt,  if  we  now  take  the 
limit  as  z  ->  Zp.  we  find  that 


lim  K  * 
z  ->  Zp 


COS  0 


n(V 

n^(z  > 


(36) 


To  be  completely  rigorous  in  this  derivation,  we  should  only  integrate  by 
parts  over  an  interval  frou  z0  to  z,  where  ze  is  sufficiently  close  to  zp  to 
guarantee  that  n'(z  )  £  0  for  z#  1  z  i  Zp.  Of  course,  the  saae  result 
obtains  because  the  remaining  integral  on  the  interval  form  to  zv  remains 
finite  as  z  ->  Zp  while  the  square  root  multiplier  approaches  zero  in  that 
limit. 

The  point  zp  is  a  turning  point  for  the  ray  labeled  by  f).  On  the 
continuation  of  the  ray,  z  decreases.  To  obtain  the  extension  to  this  ray 
continuation,  we  need  only  replace  the  integral  over  the  interval  ($,»z)  by 
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a  pair  of  integrals  over  the  intervals  (z,zp)  and  (ti«zp': 


(37) 


The  determination  of  K  is  not  so  straightforward.  We  cannot  simply 


differentiate  the  first  line,  here,  to  determine  dxx/d p,  because  the 
resulting  integral  diverges.  Ve  must  first  integrate  by  parts  as  was  done 
above  to  analyze  E  near  z 


and  then  differentiate  with  respect  to  fi.  This  leads  to  the  following 
result  for  I: 


I  =  l|n*(z)  -  sin*  fi 

n*^)  -  2  sin*  fi 
nUt)n'U,)  ^n*({t)  -  sin*  fi 


n*(z)  -  2  sin*  g 


n( 


z)n'(z)  ^n*(z)  -  sin*  fi 


ft.  *  1 

t 

N 

+ 

N  * 

XD 

1  ]  n*(z ')  -  2  sin*  0 

ln(z ')n'(z') J  ^n*(x')  -  sin*  fi 

(39) 


Again,  we  Bust  take  care  not  to  use  this  result  if  n'(z)  «  0  over  some  part 
of  the  doaain.  Ve  would  revise  this  result  by  using  the  integration  by 
parts  technique  only  over  a  portion  of  the  doaain  of  integration  near  z^  in 
(37)  and  siaply  differentiating  with  respect  to  fi  in  the  reaaining  integral. 

We  reaark  that  now  1  can  be  zero  for  an  appropriate  c(z).  That  is, 
beyond  the  turning  point,  the  ray  faaily  can  have  a  caustic. 

This  coapletes  our  discussion  of  the  in-plane  analysis  of  the  21/ ,d 
asymptotic  Green's  function. 


4.  H«  Ir-Hut  Kirehhoff  ^pioiiutin  is  2V,D. 


Ve  shall  now  consider  the  following  problem  in  s  c(x,  z)  medium.  A 
point  source  is  located  at  (tx>0,  (,).  a  receiver  is  located  at  an<* 

a  reflecting  surface*  S,  is  located  at  depth.  Consistent  with  onr  interest 
in  21/aD,  we  assme  that  the  reflector  is  cylindrical/  that  is.  the 
reflector  is  generated  by  a  curve  in  the  (x.  z) -plane  and  straight  lines 
through  that  curve  parallel  to  the  y-axis.  Thus,  we  characterise  the 
reflector  by  the  generating  curve  C  in  the  (x, z) -plane. 


C:  x  -  x^ts).  z  *  x^s) 


Ve  shall  take  s  to  be  an  arclength  variable  along  C.  Our  objective  is  to 
obtain  an  asymptotic  representation  of  the  upward  propagating  wave  arising 
in  response  to  the  reflector  under  the  assumption  that  the  downward 
propagating  wave  has  no  caustics  before  impact  with  the  reflector.  Our 
representation  is  to  be  given  totally  in  terms  of  in-plane  functions. 


The  reflector  is  to  be  a  surface  of  discontinuity  of  the  velocity 
profile,  with  no  discontinuities  above  it.  Thus,  we  will  denote  by  c_(s) 
end  c+(s)  the  limits  of  c(x, z)  from  sbove  and  below  S.  respectively. 


Ve  begin  our  analysis  with  the  Kirehhoff  integrsl  representstion  of  the 
upwsrd  scattered  field  [MMVP(8.4.S>] : 
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If  ao(*/a)  «»-<*'£>  1 

[  V5*’^  sr—  -  G(S'*a)  “fs - J  «  •  <«> 

s 

In  this  equation,  Ug(gf£)  denotes  the  response  at  a  due  to  a  sonrce  at  £/ 

G(z;g)  is  the  free  space  Green's  function  at  z  dne  to  a  sonrce  at  g*  if  Bn  is 

the  upward  normal  derivative  on  S. 

Ve  shall  stake  a  number  of  asymptotic  approzimations  in  the  inte¬ 
gral  (41).  First,  we  shall  use  for  G  the  asymptotic  Green's  function  of 
Section  2.  We  cannot,  as  yet,  use  the  results  of  Section  3,  since  the 
integration  in  (41)  is  over  a  domain  which  eztends  ont  of  the  plane,  y  *  0. 
To  emphasize  the  dependence  on  the  source  point,  we  write  the  solution  (2) 
as 


G(zsg)  -  A(g/g)  ezp  {iwc(g/g))  .  (42) 

For  the  normal  derivative  of  G,  we  only  calculate  the  leading  order  term  at 
high  frequency,  arising  from  differentiating  only  the  ezponent: 

— g— - i«l,VT(z/g)  M  z/g)  ezp(iwr(z/g) )  (43) 


Here,  £  is  an  upward  unit  normal  vector  and  the  gradient  is  to  be  calculated 
with  respect  to  the  variables,  g.  It  would  be  inconsistent  to  retain  the 
term  arising  from  the  differentiation  of  A,  itself,  since  this  is  of  the 
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tia*  (lover)  order  es  e  term  arising  froa  a  first  order  correction  to  A 
Multiplied  by  i<*fi«Yx. 


For  ug(x;£)  on  S.  ve  will  nse  a  ray  theoretic  approximation.  This 
field  is  the  reflected  wave  on  S  due  to  the  incidence  of  the  wave  in 
response  to  the  point  source  at  £.  say.  uj  (x;£).  This  function  is  also  a 
Green's  function: 


*!<?'£>  “  G(x/£) 


The  value  of  Ug(xs&)  on  S  is  just  the  ray  data  for  the  upward  reflected 
wave.  This  ray  data  is  derived  in  [)MWP,  Section  8.31.  We  find  that 


us(x.*£)  -  Ruj(x;£)  ,  x  on  S  , 


where  the  ray  theoretic  reflection  coefficient  is  given  by  [MNWP(8.3.47>] 


dxT/dn  -  Isgn  dxT/dn|  Jc  *-  c  +  ldxT/dn] 


axj/an  +  [.gn  axj/dnj  |  c-*_  c-»+I,tiy 


In  this  equation,  we  have  used  the  notation 


tfTI 

tj  *  x(x;£)  /  yj—  ■  S’VxCi/i) 


Og(i;£)  ~  R  A (*;£)  exp  {i*rc(j;£))  • 

dn«(x;£.) 

- - - iwft*VT(x;£)R  A(j;l)  exp  [iar(s;£))  . 


(48) 


The  ainus  sign  in  the  second  line  here  is  not  obvious.  It  crises  froa  tike 
feet  that  the  norasl  coaponent  of  the  gradient  of  the  phase  on  the  reflected 
wave  must  be  opposite  in  sign  to  the  noraal  coaponent  on  the  incident  wave 
in  order  that  the  incident  and  reflected  waves  be  oppositely  directed  with 
respect  to  S. 

These  two  approxiaations  constitute  a  general ixat ion  of  the 
approx iaat ions  credited  to  Rayleigh.  Soaaerfeld.  Fresnel  and  Rnygens.  as 
well  as  Kirchhoff.  For  discussions  of  the  distinctions  between  these 
approxiaations.  see  Ooodaan  [1968].  Kuhn  and  Alhilali  [1977]  and  Volf  and 
Marchand  [1964]  .  For  siaplicity.  we  shall  refer  to  (46)  as  the  Kirchhoff 
spproxiaation. 

Ve  substitute  (40),  (41)  and  (46)  into  (31)  and  obtain 


Ujg^.'a)  -  i«  J  RA(g;£)  A (j/n)  exp{ia[t(j;J)+r(x;a)]} 
S 

•  i‘[V-r(g/T|)  +  Tt(x/£>]  dS 


(49) 


This  asyaptotic  rsprssentatioa  Is  givea  la  tens  of  the  2*/*D  ray  theoretic 
aaplitade  aad  phase  iatrodaced  ia  Sectioa  2.  However,  it  still  requires  oat 
of  plaae  coapatatioas  of  aaplitade  aad  phase. 


Ve  will  aow  reduce  the  represeatatioa  (49)  to  oae  ia  which  the 
iategraad  is  givea  totally  ia  tens  of  the  ia-phase  aaplitade  sad  phase.  To 
do  this,  we  eote  first  that 


dS  *  dsdy  =  ds  dij  .  (SO) 

That  is.  the  differeatial  surface  area  is  the  product  of  the  differeatial 
arc  leagth  ds  aloag  the  geaeratiag  carve.  (40),  aad  the  differeatial  out-of¬ 
phase  coordiaate. 

We  will  calculate  the  iategral  ia  y  «  z,  by  the  aethod  of  statioaary 
phase  [MMWP,  Chap.  2] .  The  ]Aase  is 


f( J/ J/3)  =  t(s;£)  +  ^(5/3)  » 


(51) 


with  first  derivative 


el 


ar(j.'i)  aT(i/n) 


P,U.’£)  +  P1<5'3) 


(52) 


la  the  last  eqaatioa,  ve  have  used  the  aotatioa  of  (7),  bat  with  arguaeats 
that  properly  relate  each  p,  to  the  correspoadiag  x. 


At  the  stationery  point,  di/9xa  -  0.  Thus,  the  an  of  the  pa’s  neat  be 
zero.  Froa  (10)  and  (11)  ve  see  that  both  pa'a  will  have  the  saae  sign  as 
i,,  for  xa  A  0.  Thus,  the  sna  can  only  be  zero  when  x,  ■  0.  That  ie,  the 


phase  has  a  stationary  point  in  xa  only  for  xa  “  0.  In  tema  of  the 
paraaeters  a  and  0,  this  oeenrs  at  a  *  0. 

Ve  anst  now  determine 


a*f 


9x. 


«pa<i 3pa<g/a) 

-  +  - 

dx.  dx. 


(53) 


it  i,  *  0  or  «  *  0.  To  do  this,  we  rewrite  (11)  as 


xs  “  P.  * 

and  differentiate  iaplicitly  with  respect  to  x„: 


(54) 


1 


(55) 


e  at  pa  * 


When 


lnat 


0,  we  find  that 


3x,  a  *  P* 


Tki«  ratal t  applies  equally  to  p,(*/£)  and  p»(*/n),  except  that  o  aust  be 
properly  interpreted  as  tbe  paraaeter  on  the  corresponding  ray/  that  is,  the 
paraaeter  connecting  x  to  £  or  x  to  a.  To  distinguish  between  these  two  ray 
paraaeters,  we  shall  denote  tbe  first  by  eg  and  the  second  by  Thus, 


J_  +  JL 


Ve  note  that  sgn  df/dxa  »  +  1  at  x,  =  0.  This  result  is  needed  for  the 
stationary  phase  evaluation.  Furthermore,  the  entire  integrand  is  to  be 
evaluated  at  xt  “  0,  that  is,  in-plane .  Thus,  we  nay  use  the  results  of 
Section  3,  to  determine  A,  t  and  o,  now  for  each  of  the  fixed  points,  J,  and 
ij.  The  result  is 


ug  (£,*a>~ 


|2w|wl  exp{3wi/4  sgnu} 


R  exp  (i«CT(x,'i)+T(x,*a^ I 


C,xa=  0 


hwr*  *nX) 


•  B*[Vt(x;£)  +  VT(x,*a>]  ds. 
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N  %  , 


Ia  this  equation,  if*  and  are  the  in-plane  Jacobiana  defined  by  (17)  for 
the  rays  eaanating  fro*  J  and  3.  respectively.  Evaluation  of  the  integrand 
»t  x,  *  0  reduces  all  of  its  elesients  to  the  in-plane  results  of  Section  3. 

It  is  also  worthwhile  to  specialize  (58)  to  the  case  of  backscatter  or 
zero-offset,  ug(£;£).  We  can  then  diapense  with  the  subscripts  (  and  q. 
The  result  is 


* 

_  (i/S)  -  J*g»wi/4  sgnw)  Reap  (2i«T(x/&)) 

8n  '*  K  fT 


£)  ds 


C.  zf-  0 


(59) 


The  results.  (58)  and  (59)  are  asyaptotic  representationa  of  the  21/SD 
in-plane  reflection  response  for  non-zero  or  zero  offset  observations, 
respectively.  These  foraulas  are  starting  points  for  either  nuaerical 
coaputation  of  the  reaaining  integral  or  further  asyaptotic  analysis  —  say, 
stationary  phase  in  the  arc  length  s.  Thus,  we  have  achieved  a 
representation  of  the  in-plane  fields  in  the  2*/aD  case  with  the  properties 
carried  out  in-plane  (xa  »  0)  and,  nonetheless,  the  effect  of  geoaetrical 
spreading  is  accounted  for,  asyaptotically. 

As  in  Section  3,  we  now  consider  the  special  cases  of  constant 


background  velocity  and  depth  dependent  background  velocity.  For  the  foraer 
case,  we  use  the  results  (23,  (24)  and  (25)  to  obtain  in  place  of  (58) 


frj 


«s<5*a>  ~ 


hnlttl  exp  {-nil 4  sgn  «} 


(4n)*  ^jc# 


R  exp  (i«[r  +  r  ]/c  } 


r  r  l-i  -i 

**  *r£  +  % 


C.x%*  0 


n*[r.  +  r  ]  d» 
{  W 


la  this  equation,  rg  is  the  distance  defined  by  (22)  while  rtj  is  the  same 
function  with  £  replaced  by  3.  Furthermore,  8  is  the  unit  upward  normal  to 


i,  «  (£4-  xx,£f-  xf)/r.  =  (1^-  x^  q#-  x^/r  =  -  coVr(x.-i})  (61) 


and,  from  (46). 


a-*  /co  -  d.+-  -  <r 


T*  -  c~*  +  [n'?./c  ]* 


+  L'  -  C  1  +  [&*ipfc  ] 
5  0  1  +  0  5  0 


c  =  c 
0  - 


We  note,  here,  that  8  and  fg  are  colinear  and  sgn  8*fg  =  +1.  It  is  apparent 
from  this  resnlt  that  we  have  momentarily  lost  the  symmetry  of  ug  in  £  and 
3.  However,  were  we  to  apply  the  method  of  stationary  phase  to  (60),  we 
would  find  that  the  dominant  contributions  to  the  field  at  £  arise  from  the 
specular  points,  at  which  8*8^  «  8’8g  and  the  symmetry  would  be  restored. 
Alternatively,  we  could  redefine  (62)  as  a  reflection  coefficient  Rg  to 


-.1 

*  .  •  .  1 


,*.y  •  yV  .'A*.  *  V 

*,*  *  "  *,*  »,*  *,*  *  *  o  -  *  *,■  •  *  • 


.  '  »  *  1  « 

,  h  ‘  * i 


with  r  given  by  (23) 


Ve  now  turn  to  the  cate  of  a  depth  dependent  background  velocity.  In 
this  case,  we  can  do  little  more  than  substitute  into  (58)  or  (60)  the 
parametric  representations  given  by  (32)  and  (34)  —  one  for  the  ray  family 
emanating  from  £  with  ray  parameter  another  for  the  ray  family  emanating 
from  tj  with  ray  parameter  0_. 


This  completes  the  discussion  of  the  asymptotic  2*7,0  in-plane  upward 


scattered  field. 


5.  H<  Affect  of  Variable  Density 


It  is  fairly  straightforward  to  include  the  affects  of  variable  density 
in  the  results  of  the  previous  sections.  Variable  density  does  not  change 
the  phase  of  the  solntion.  bnt  only  the  aaplitnde.  Following  Brekhovshikh 
[1980],  if  we  assume  that  n  denotes  pressure,  then  the  governing  equation 
(1),  is  replaced  by  (op.  cit.,  eq.  (19.2)) 

* 

pV-  [(l/p)Vu]  +-!^u  =  -  V  &(*,)  S(x,-  $,)  .  (64) 

c 

In  this  equation,  p  *  p(x)  is  the  variable  density. 

Ve  will  again  assume  that  u  is  of  the  form  (2).  However,  in 
anticipation  of  a  change  in  amplitude  (but  not  in  phase)  we  shall  denote  the 
amplitude  by  B: 


u(x/w)  ~  B(x)  exp  {i«  t(x) }  .  (65) 

The  phase  r  is  again  a  solution  of  the  eikonal  equation  (3),  but  B  is  a 
solution  of  the  transport  equation 

2VB-Vt  +  BpV( 1/p) *Vt+  BV*t  =  0  .  (66) 

If  we  divide  here  by  p1^*,  we  find  that 
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2V(Bp  '  ) *Vt  +  Bp  '  Ti=0  .  (67) 

That  is.  Bp-1/1  satisfies  the  same  transport  eqoation  as  A  does  —  equation 
(4)  —  and  hence. 


B(x)  =  const.  A(x)  ip(x)  .  (68) 

Here,  const,  means  constant  with  respect  to  the  running  parameter  along  the 
rays.  In  fact,  if  one  were  given  the  same  ray  data  for  B  as  we  were 
previously  given  for  A,  then  the  initial  value  of  B  would  have  to  agree  with 
the  previous  initial  value  of  A.  From  this,  it  follows  that 


B(x)  -  A(x)  Jp(x)/p(x#) 
where  x,  is  the  initial  point  on  the  ray. 


(69) 


The  Green's  function,  developed  in  Section  2  is  precisely  a  solution 
for  which  the  initial  data  for  A  at  x  =  J  —  equation  (6)  —  now  becomes  the 
initial  data  for  B.  Thus,  from  (12)  we  immediately  determine  that  for  3D 
propagation 


Jp(x(o,  u,  P>)  sin 
B  - - 1  ■  ■  — 

4jt  1c($)  p(£)  J(e, a, 0) 

with  J  again  defined  by  (13). 


(70) 
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la  the  specialization  of  Section  3  to  the  in-plane  21/aD  Green's 
function,  we  find  froa  (18)  that 


j  p(z(a, 0)) 

4jt  '  a  Zp(£^ 


(71) 


with  K  again  given  by  (17)  and  z(n#P)  the  solution  of  the  ray  equations  in 
(14).  Purtheraiore.  the  representation  (19)  is  now  replaced  by 


1  |p(z(o.  0)) 

n  ~  4?  fTYTTIT  exp  lim]  *  (72) 

The  amplitude  B  can  be  related  to  the  amplitude  Ba  0f  the  two- 
dimensional  asymptotic  Green's  function  by  (21)  with  A's  replaced  by  B's. 
However,  Ba  is  related  to  Aa  by  (69)  with  subscripts  2  introduced  on  both 
sides  of  the  equations. 

For  the  case  in  which  p  and  c  depend  only  on  z  or  z(>  we  use  (32)  for  t 
and  (34)  or  (35)  for  K. 

These  results  are  to  be  substituted  into  (72)  along  with 


p(z(o, 0) ) 


p(z) 

pTTJ  * 


(73) 


The  discussion  of  the  in-plane  Kirchhoff  approximation  in  21/,D, 


presented  in  Section  4,  anst  also  be  revised  to  allow  for  a  jump  in  density 
as  well  as  a  jnap  in  velocity.  The  reflection  coefficient  given  in  (46)  was 
derived  nnder  the  assumptions  that  n  and  du/dn  are  continnons  across  S.  For 
the  differential  operator  (64),  we  anst  replace  those  conditions  by  the 
conditions  that  n  and  (1/p)  dn/dn  anst  be  continnons  across  S.  Proceeding 
as  in  the  derivation  of  (46)  given  in  [MNWP,  Section  8.3],  we  deduce  that 
the  proper  reflection  coefficient  is  now  given  by 


p_*  dtj/dn  -  p+*  (sgn  dx^/dn]  ^c+*  -  c_*+  [dXj/dn]' 

R  - -  .  (74) 


p1  dTj/dn  +  p+J[sgn  dx  £dn]  Jc+*  ~  c_ 


+*  -  c_*+  [dXj/dn]* 


In  (49),  we  anst  replace  A(x;£)  A(z;i))  by  B(x,*£)  B(x/r)  and  nse  (74) 
for  R.  With  these  adjnstaents  (38)  is  replaced  by 


For  the  cue  of  i  depth  dependent  density  function,  we  need  only 
replace  p(£)»  p(g)  and  p(z)  by  p((j)  and  p(gt)  and  p(zt),  respectively. 
Furthermore,  we  nse  (32)  and  (34)  for  t,  o,  and  K  while  still  using  (74)  for 
R.  These  results  can  be  nsed  in  (7S)  or  (76). 

This  completes  our  discussion  of  the  Variable  Density  case. 


6 .  CONCLUSIONS 


A  fev  results  aboot  wave  propagation  in  two-and-one-half  dimensions 
have  been  derived.  The  first  set  have  to  do  with  the  asymptotic  radiation 
of  acoustic  waves  from  a  point  source.  The  second  set  are  related  to  the 
Kirchhoff-approximate  upward  scattered  field  from  a  cylindrical  surface  at 
depth.  In  both  cases,  we  have  found  that  the  in-plane  propagation  of  a  wave 
in  three  dimensions  can  be  described  totally  in  terms  of  in-plane 
calculations  which  are  no  more  difficult  than  would  be  carried  out  to 
generate  two  dimensional  models. 

In  our  inversion  research,  we  regularly  develop  our  theory  as  if  the 
observations  are  known  over  a  planar  array  on  the  upper  surface.  However, 
we  then  specialize  our  results  to  the  2‘/,D  case  in  which  data  is  known  only 
on  a  line.  Ve  have  already  begun  using  the  results  described  here  to 
generate  synthetic  data  to  test  these  2*/,D  algorithms.  This  was  the 
motivation  for  this  project  and  the  utility  to  our  group  has  already  been 
demonstrated. 
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TABLE  I 


* 


c(x, z) 


19. 17.14 


IN-PLANE  GREEN'S  FUNCTION 


CONSTANT  DENSITY 


c(z) 


19.32.34 


2. 24,  26 


VARIABLE  DENSITY 


72,17.14 


72,32,34.73 


t 


TABLE  II 


nonzero 

offset 

zero 

offset 


nonzero 

offset 

zero 

offset 


t  *  v*v 


IN-PLANE  KIRCBHOFF  APPROXIMATE  FIELD 


CONSTANT  DENSITY 


c(x, z) 


c(z) 


c 


0 


58,46.17.14  58.46,34.32  60.61.62,23 


59.46.17.14  59,46.34,32  63.61,62,23 


VARIABLE  DENSITY 


75.74. 17,14 


75,74.34,32 


76,74.17.14  76,74,34,32 
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